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Abstract 

We propose a general class of five-dimensional soft-wall models with AdS 
metric near the ultraviolet brane and four-dimensional Poincare invariance, 
where the infrared scale is determined dynamically. A large UV/IR hi- 
erarchy can be generated without any fine-tuning, thus solving the elec- 
troweak/Planck scale hierarchy problem. Generically, the spectrum of fluc- 
tuations is discrete with a level spacing (mass gap) provided by the inverse 
length of the wall, similar to RSI models with Standard Model fields propa- 
gating in the bulk. Moreover two particularly interesting cases arise. They 
can describe: (a) a theory with a continuous spectrum above the mass gap 
which can model unparticles corresponding to operators of a CFT where 
the conformal symmetry is broken by a mass gap, and; (b) a theory with 
a discrete spectrum provided by linear Regge trajectories as in AdS/QCD 
models. 



1 Introduction 



Warped extra dimensions were introduced by Randall and Sundrum [1] (RSI) as 
a very elegant way of solving the hierarchy problem by means of the geometry of 
extra-dimensional theories. The original proposal consisted of a slice of AdS space 
bounded by two branes: the ultraviolet (UV) brane, located closely to the AdS 
boundary, and the infrared (IR) one. One of the most exciting aspects of RSI 
theories is provided by the AdS / CFT correspondence [2] , by which fields living in 
the UV brane are fundamental fields which can interact with a strongly coupled 
conformal sector (CFT), while fields living in the IR brane are the dual description 
of operators of the CFT. In this way a Higgs boson localized on the IR brane 
provides a dual description of the composite Higgs as a bound state of an extra 
strong interaction (technicolor). From the five dimensional (5D) point of view the 
Higgs boson mass is simply redshifted from its natural value, of the order of the 
Planck mass, to the TeV scale by the warp factor. The IR brane spontaneously 
breaks the conformal symmetry and its location needs to be stabilized by some 
dynamical mechanism. 

A dynamical mechanism to stabilize the IR brane was proposed by Goldberger 
and Wise [3] (GW) by the introduction of a background scalar field propagating 
in the 5D bulk and acquiring a coordinate dependent vacuum expectation value 
(VEV). This triggered a 4D effective potential for the radion field with a minimum, 
stabilizing the brane separation at scales of order 1/TeV, determined by the values 
of the scalar field at both the UV and IR branes, with a modest ~ 1% fine- 
tuning. Furthermore the back-reaction of the scalar field through the 5D Einstein 
equations generated a deviation from AdS of the metric far away from the UV 
brane, still preserving the main features of the AdS / CFT correspondence. In fact 
this is a general feature since AdS is the only 5D metric consistent with no (or 
constant) background scalar field, which means that any stabilizing bulk field is 
expected to back-react on the AdS metric. 

Since the presence of bulk scalar fields, required to stabilize the brane dis- 
tance, is expected to disturb the AdS metric in the far IR, one can introduce 
phenomenological models with background scalar fields and with the only con- 
straint of describing the AdS geometry near the UV brane. One of them is of 
course RSI with GW mechanism, but there are more general models where the 
IR brane is replaced by a naked curvature singularity [5-8]: these models are 
called soft-wall models [9-14]. Soft- wall models possess AdS geometry near the 
UV brane and no IR brane, i.e. they have a non-compact extra dimension as those 
proposed by Randall and Sundrum with an infinite extra dimensional length [15] 
(RS2), but with additional background scalar fields. These fields back-react on 

1 For an earlier analysis of models with exponential (Liouville) dilatonic brane potentials see 
Ref. [4]. 
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the metric and generate a singularity at a finite value of the extra dimensions: 
the (finite) length of the extra dimension provides the location of the soft wall. 

Soft wall models were first introduced to model the Regge behaviour of ex- 
cited mesons in AdS/QCD models [9,10], as an alternative to RSI for electroweak 
breaking models solving the hierarchy problem and providing experimental signa- 
tures at LHC [12,13] and also to provide a 5D setup to the theories of unparticles 
recently proposed by Georgi [16] in the presence of a mass gap [14, 17]. In all 
these theories, as in GW models, there are background bulk scalar fields and thus 
a built-in mechanism to stabilize the extra dimension, i.e. the distance between 
the UV brane and the singularity. However in the proposed models this distance 
is naturally of the order of the AdS length and the hierarchy problem does not 
find a natural solution, i.e. for values of the scalar field at the UV brane of the 
order of the 5D Planck scale [13]. 

In this paper [18] we will propose a set of one-parameter (u) 5D models with a 
background scalar field propagating in the bulk of the extra dimension and with 
the following properties: 

• The metric is AdS near the UV brane. 

• The hierarchy between the AdS length and the soft-wall length (the elec- 
troweak length) is naturally stabilized for values of the scalar field at the 
UV brane of the order of the 5D Planck scale. The naturalness of this result 
comes from a double exponential suppression. 

• For values of v < 2 the naked singularity does not contribute to the vacuum 
energy and thus it does not need to be resolved to satisfy Einstein equa- 
tions. Moreover the potential is bounded from above in the solution and the 
singularity is a physical one according to Ref. [7] . 

• For 1 < v < 2 the spectrum of fluctuations behave as in RSI models with 
Standard Model (SM) fields propagating in the bulk: for scalar field fluctu- 
ations localized near the singularity the hierarchy problem is automatically 
solved. 

• For v = 1 one can model unparticles with a mass gap provided by the inverse 
length of the extra dimension. Fluctuations have a continuous spectrum 
above the mass gap. 

• For < v < 1 unparticles without a mass gap emerge: the spectrum is 
continuous above zero mass. 

The paper is organized as follows: In Section [2] we review the general construc- 
tion of backgrounds and the associated question of self vs. fine-tuning of the cos- 
mological constant. The general conditions for models with physical singularities 
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are summarized. In Section [3] we consider the background solutions for a particu- 
larly simple class of models satisfying all physical requirements on the singularity 
and describing the electroweak/Planck hierarchy without fine-tuned parameters. 
Fluctuations on the background for the graviton, radion and the scalar field are 
studied in Section H] where a complete numerical analysis and some analytical ap- 
proximations are provided. A general class of models having all the good required 
properties, included the hierarchy determination, is presented in Section [51 One 
of these models has the mass of excitations as m 2 ~ n and it is thus a good candi- 
date to model the Regge behaviour in AdS/QCD models. Finally our conclusions 
and outlook are drawn in Section [6j 



2 The 5D Scalar Gravity System 

In this section we would like to review the construction of backgrounds with 4D 
Poincare invariance proposed in Ref. [19,20], and the associated question of self 
vs. fine-tuning of the cosmological constant (CC) [5-8,21,22]. Let us consider 5D 
gravity with a scalar field, and look for the most general solutions to this system 
that preserve 4D Poincare invariance, i.e. a background of the form 

ds 2 = e-^dx^dx^ + dy 2 , (2.1) 

with a "mostly plus" flat Mikowskian metric rj^ = diag(— ,+,+,+, +) and an 
arbitrary warp factor A(y). We will introduce a single brane sitting at y = and 
impose the orbifold Z 2 symmetry y — > —y under which A and <fi are even. We will 
thus consider the bulk plus brane action 

S = Jd 5 x^ [M 3 R - 3(«90) 2 - V ((/>)] - J d 4 x V /3 ^A(0). (2.2) 

We have introduced arbitrary bulk and brane potentials V and A, and M denotes 
the 5D Planck mass which we will set to unity for the remainder of the paper. 
Note the noncanonical form of the kinetic term for the scalar that will simplify 
future formulas. The bulk equations of motion (EOM) that follow from the action 
in Eq. (J22D are 

6(f)" -2AA'(f)' -d^V ((f)) = 0, (2.3) 
A" -(f)' 2 = 0, (2.4) 
12A' 2 - 3(f)' 2 + V = 0. (2.5) 

This system has three integration constants 0. One of them is A(0) that remains 
totally free. The other two can be fixed from the boundary conditions (BC) that 

differentiating the third equation one obtains a linear combination of the first two, such 
that (say) the first can be discarded. The system is then first order in <p and second order in A. 
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follow from the boundary pieces of the EOM, 

A'(0+) = -U(0 O ), (2.6) 
o 

0'(O + ) = ^A(0o), (2.7) 
o 

where 0o = 0(0). Using Eqs. (I2.6P and (12 .7p in Eq. (12.51) determines 0o, 

^A(0 O )] 2 - ^A(0 O ) W(0 O ) , 

which could be used to replace (12.71) . 

The authors of Ref. [20] introduced the following "trick" to obtain solutions to 
this system by defining the so-called "superpotential" via the differential equation 

3(d^W) 2 - 12W 2 = V, (2.9) 

and writing 

A' = W(<j>), 
0' = W), 

while the boundary conditions are satisfied if 

W(<h) = ^A(0 O ) , d+Wifa) = ^A(0 O ). (2.12) 
o o 

Again, the system of Eqs. (I2.9I) - (I2.11I) has three integration constants and in 
principle every solution to Eqs. (I2.3p - (I2.5I) can be constructed in this way [20]. 
One integration constant is the trivial additive constant ^4(0) that does not enter 
in Eq. (1 2 . 1 2 [) . We are left with the integration constant in Eq. (I2.9p and the value 
O to fix the two constraints Eq. (12.121) . The equation for W is a complicated 
nonlinear differential equation, and in practice it is often easier to start with 
a particular superpotential satisfying the boundary conditions and deduce the 
potential needed to reproduce it. 

A peculiarity of the scalar-gravity system with one brane is the appearance of 
naked curvature singularities at finite proper distance. In particular, it can easily 
be checked that if the superpotential W grows faster than 2 at large 0, the profile 
4>{y) diverges at finite value of y = y s . Moreover the 5D curvature scalar along 
the fifth dimension can be written as 

R(y) = 8(d^(0[y])) 2 - 2OW(0[y]) 2 , (2.13) 

so that the curvature in general diverges at y — y s . The interpretation is that 
spacetime ends at y s . 
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(2.8) 



(2.10) 
(2.H) 



Having three integration constants but only two constraints, it seems that 
one can obtain flat 4D solutions with fairly generic brane and bulk potentials 
without fine-tuning. This miraculous self-tuning property of the scalar-gravity 
system was first pointed out in [5, 6] and it was further scrutinized in several 
papers [7,8,21,22]. In particular, the authors of Refs. [21,22] pointed out that 
the on-shell Lagrangian, integrated over the fifth dimension can be written as 

£° n - shcU = \K<Po) + | Jdy e-^V^y]) (2-14) 

(we have set A(0) = 0). They then make particular choices for V and A and show 
that the result is non-vanishing. The interpretation of this apparent contradiction 
to the existence of a flat background is simple: having dynamically generated a 
new boundary at the singularity, we must ensure that the boundary pieces of 
the equations of motion vanish at y = y s . If this is not the case, the resulting 
"solution" does in fact not extremize the action, resulting in a nonzero 4D CC. 
This can actually be seen in rather general terms. Making use of the equations 
for the superpotential, we write 

e-' A ^V{<P[y]) = 3 A{e- 4A ^W/(<%])}, (2-15) 

leading to 

£ on - shc11 = ^A(0 O ) - 2W{<k) + 2 lim e-^W^y}) . (2.16) 
o y^y s 

The first two terms cancel if Eq. (12.121) is satisfied, while the last one depends on 
the particular form of the superpotential. In order for the last term to vanish, 
W needs to grow more slowly than e 2 * at large 0: this can be seen by using the 
field itself as a coordinate. The position of the singularity moves to <f)(y s ) = oo 
and the equation for A(<f)) becomes A' = W/W. It follows that the last term in 
Eq. (12.161) goes to a constant for W ~ e 2< ^, while it goes to infinity (zero) when 
W grows faster (slower) than e 2r ^. We thus arrive at a simple criterion for the 
existence of singular solutions: 

A singularity with 4>(y s ) —> oo is allowed if and only ifW{4>) / 9 17 % 

grows asymptotically more slowly than e 2 ^. 

Notice that a potential (V) growing more slowly than e 4< ^ is only necessary but 
not sufficient for the validity of (j2.17|) . the trivial counterexample being V = 
which has the general solution W = ce 2 ^. It is instructive to compare our 
criterion with the one found in Ref. [7] where AdS-CFT duality was used to 
classify physical singularities. According to Ref. [7] admissible singularities are 
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those whose potential is bounded above in the solution. Inspection of Eq. (12. 9p 
shows that singularities fulfilling (I2.17P have a potential that goes to — oo, while 
those that fail (12.171) go to +00. Although we here employ a much more basic 
condition (a consistent solution to the Einstein equations), which in particular can 
be applied to theories without any field theory dual, it is good to know that our 
allowed solutions have potentially consistent interpretations as 4D gauge theories 
at finite temperature. 

However achieving a superpotential that grows more slowly than e 2 ^ requires a 
hidden fine-tuning of the cosmological constant. To see this it suffices to consider 
a potential that behaves asymptotically as 

V~be 2 ^. (2.18) 

Writing W as 

W{(j>) = w((p)e^, (2.19) 
we can express the solutions for w as the roots of 

e (4-^)(*-c) = ( 2w + V6 + 4w 2 ) ±2 (i/w =F Vb + Aw 2 ) u , (2.20) 

where c is an integration constant. For v > 2 this implies that w asymptotes to 
a constant 

(2 - 2i) 

at large (j) an d for b > 0. However for < v < 2, w generically behaves as 

w-e^ 2 -^. (2.22) 

Only if we adjust c — > 00 we can achieve that w behaves as in Eq. (12.211) . In this 
case, b has to be negative in order to have a real solution for W. 

The generic solution to Eq. (12. 9p thus grows as W ~ e 1 "^ for v > 2 and W ~ e 2 ^ 
for v < 2. However, it is possible to arrange for W ~ e v< ^ in the latter case by 
picking a particular value for the integration constant in Eq. (12.91) [§. There are 
thus two possible scenarios. 

• The superpotential W grows as e 2 ^ or faster and the equation of motion are 
not satisfied at the singularity. The only consistent way out is to resolve 
the singularity, for instance by introducing a second brane located at y s . 
In that case fine-tuning of the CC is restored as we introduce two more 
conditions analogous to Eqs. (12.61) and (12.71) or Eq. (12.121) respectively, but 
do not increase the number of free parameters [21]. 

3 Similar reasonings apply to potentials that grow even more slowly, e.g., as a power. The 
generic solution behaves as e 2< ^, but particular solutions may exist that behave as VV and hence 
allow for consistent, yet fine-tuned, flat backgrounds. 
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• The superpotential grows as e u ^ with v < 2, or slower. The equations 
of motion are satisfied at the singularity, and there is no need to resolve it. 
The price one pays is the adjustment of the integration constant in Eq. (12.91) . 
In this case we lose one of our parameters needed to satisfy the boundary 
condition Eq. (12.121) . resulting in a fine-tuning of the brane tension. 

It is important to realize that either fine-tuning precisely corresponds to the fine- 
tuning of the CC. In the second possibility above this is particularly obvious: the 
superpotential is completely specified by the bulk potential and the boundary 
condition at — > oo. Eq. ( 12.121) is then simply the minimization of the 4D 
potential 

V A ((f>) = \{<f>) - 6W(<p) (2.23) 

under the condition that V^(0) vanishes at the minimum = O - in fact the 
brane potential A(0) should be determined by physics localized at the UV brane 
interacting with the (dilaton) field <fi. For example if the Standard Model Higgs 
field H is localized at the UV brane it will generate a brane potential as \(<f), H) 
which will in turn provide the effective brane potential X(4>) = A(0, (H)) after 
electroweak symmetry breaking. So after the electroweak phase transition □ there 
will be a ^-dependent vacuum energy which will require re-tuning the cosmological 
constant to zero and possibly a shift in the minimum of Eq. f l2.23[) . 

What matters to us here is that there exist consistent solutions to the equations 
of motion in the full closed interval [0,y s ] that, although requiring a fine-tuning 
of the CC, do not demand the introduction of a second brane or any other means 
of resolving the singularity. 

3 A Model of Soft- Wall Stabilization 

Our goal in this section is to find solutions to the 5D scalar-gravity system that: 

1. require only one brane at y — 0, 

2. behave as AdS§ near that brane, 

3. have a soft wall at a value y = y s , i.e. give rise to finite volume, 

4. possess a mass gap or level spacing hierarchically smaller than the Planck 
scale without fine-tuning of parameters. 

The last requirement might be called the problem of soft-wall stabilization. We 
will now show that all these requirements can be met by starting from the simple 
superpotential 

W = k(l + e u *), (3.1) 
4 Othcr phase transitions, as e.g. the QCD phase transition, should have a similar effect. 
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where k is some arbitrary dimensionful constant of the order of the 5D Planck 
scale, and v < 2. The solution can be immediately written down 

A(y) = %-llogfl-^V 

v 2 \ VsJ 
<l>(y) = -Uog[v 2 k{y s -y)}. (3.2) 

At the point y = y s we encounter a naked curvature singularity as explained in 
Section [2j For y <C y s , i.e. near the boundary at y — 0, the geometry is AdS$. 
The bulk potential which corresponds to the superpotential (13.1ft is given by 

V(</>) = (3A;V - \2k 2 )e 2v * - 24k 2 e u * - 12k 2 . (3.3) 

• For v < 2 the potential is bounded from above. More precisely it satisfies 
the condition 

V{(j>[y\) < V{<h), (3-4) 

necessary [7] for the corresponding bulk geometry to support finite temper- 
ature in the form of black hole horizons □■ Moreover for v < 2, as we have 
seen in the previous section, the equations of motion are satisfied at the 
singularity and there is no need to resolve it. 

• For v > 2 the equations of motion are not satisfied at the singularity and the 
latter would need to be resolved to fine-tune to zero the four-dimensional 
cosmological constant. Finally the potential is not bounded from above and 
finite temperature is not supported in the dual theory. 

The location of the singularity depends exponentially on the brane value of 4>, 

ky s = - 2 e^\ (3.5) 

As we will see in the next section the relevant mass scale for the 4D spectrum is 
not the inverse volume but rather the "warped down" quantity 

p = k(ky s )~ 1/u2 e- ky ° . (3.6) 

All we need in order to create the electroweak hierarchy is thus <fto < but 
otherwise of order unity. This can be achieved with a fairly generic brane potential, 
for instance by chosing a suitable A(0) such that the second of Eq. (12. 12}) is satisfied 

5 A characteristic feature of physical singularities [7]. 
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for our superpotential 0. For negative 4>q, the ratio of scales kj p exhibits a double 
exponential behaviour 

log — + . . . , (3.7) 

p v l 

and we can create a huge hierarchy with very little fine-tuning. In Fig. [1] we plot 
p/k as a function of |0o| for different values of v and also as a function of v for 
a fixed value ky s = 30 which generates a hierarchy of about fourteen orders of 
magnitude. 

A comment about the choice of our superpotential is in order here. Its par- 
ticular form, Eq. (13.11) . guarantees full analytic control over our solution. A more 
detailed analysis of other possibilities will be postponed to Section [5j 




Figure 1: Plot of log 10 (p/k) as a function of \<f>o\ for v — 1 and v — 2 [left panel], 
and as a function of v for ky s = 30 (this value will be used in following plots) 
[right panel]. 



It will be useful in the following to define the metric also in conformally flat 
coordinates defined by the line element 

ds 2 = e- 2A ^{dx'"dx v T]^ + dz 2 ). (3.8) 

6 In order to satisfy the first of Eq. (|2.12[) we still need a fine-tuning, for instance by adding 
a <fr independent term to X(<f>). This is precisely the tuning of the 4D CC discussed above that 
of course has nothing to do with the electroweak hierarchy we want to explain here. 
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where A(z) = the relationship between z and y coordinates being given 

by exp[A(y)]dy = dz. One easily finds, for v > 0, that 

p(z - zo) = r(l - l/^ 2 , A;y s - fcy) - T(l - l/u\ ky a ), (3.9) 

where zq corresponds to the location of the UV brane that we assume to be at 
zq = 1/k and T(a,x) is the incomplete gamma function. Since we are taking 
e kys ^> 1 and hence k/p > 1 we can approximate T(l — l/u 2 ,ky s ) ~ p/k and 
(13. 9p simplifies to 

pz~T(l-l/v 2 ,ky s -ky). (3.10) 
For v > 1 the singularity at y s translates into a singularity at z s given by 

pZs ~ r(i - i/u 2 ) . (3.ii) 

For < v < 1 the singularity at y s translates into a singularity at z s — > oo. 

As we will see in the next section the case < v < 1 provides continuous 
spectra without any mass gap 0, i.e. typically it leads from a 4D perspective to 
unparticles [16]. The case v = 1 corresponds to continuous spectra with a mass 
gap provided by p in Eq. (13.71) leading in 4D to unparticles with mass gaps [14,17]. 
Finally the case 1 < v < 2 corresponds to discrete spectra with typical level 
spacing controlled by p, as in AdS models with two branes [1]. 

4 Fluctuations and the 4D Spectrum 

In this section we study the fluctuations of the metric and scalar around the classi- 
cal background solutions. A general ansatz to describe all gravitational excitations 
of the model is, with the appropiate gauge choice [23] 

<f>(? t y) = <j>(y) + <p(x,y), (4.1) 
ds 2 = e- 2A ^- 2F ^\^ v + hJJ)dxV + (1 + G(x, y)) 2 dy 2 , (4.2) 

where <fi{y) is the background solution given in Eq. (13.21) . The Einstein equations 
that arise from this ansatz have the spin-two fluctuations decoupled from the 
spin-zero fluctuations, so we can proceed to study them independently. 

4.1 The Graviton 

Let us first consider the graviton as the transverse traceless fluctuations of the 
metric 

ds 2 = e- 2A ^( VfMV + h^( Xj y))dx 2 + dy 2 , (4.3) 
7 The case v = is just the RS2 model [15] with a constant dilaton <f>. 
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where h ^ = d^W = 0. In order to respect the orbifold symmetry and to keep the 
possibility of a constant profile zero mode, we will consider h(y) = h(—y) which 
leads to the boundary condition at the brane h'(0) = 0. The part of the action 
quadratic in the graviton fluctuations becomes 

S = d^x dy \J—gR 



~+-\J d 4 xdye~ 2A(y) {d p h^dW + e- 2A{y) d y h^d y h^) . (4.4) 
Using the ansatz 

h^(x,y) = h^(x)h(y), (4.5) 

one can obtain the equation of motion for the wavefunctions h(y), which is given 
by 

h"(y) - AA'{y)h\y) + e 2A{y) m 2 h{y) = 0. (4.6) 

After an integration by parts in (14.41) . one finds an additional equation due to 
boundary terms at y = y s , 

e- 4A ^h'(y s ) = 0. (4.7) 
In addition, one has to impose that the solutions are normalizable, i.e. 



dye~ 2A ^h 2 {y) < oo. 



(4.8) 



It is now convenient to change to conformally flat coordinates, as defined 
in (13. 9p . In this frame, rescaling the field by h(z) = e~ 3A ^ z ^ 2 h(z) , Eq. (14. 6 p can 
be written clS cL Schroedinger-like equation, 

- h(z) + V h (z)h(z) = m 2 h(z), (4.9) 
where a dot denotes derivation with respect to z, and the potential is given by 

(4.10) 

The boundary equations are written in the z-frame as 



V h (z) = \A\z) - lA{z) 



-3A(z)i 



h(z) 



3A ^/ 2 (h{z) + -A{z)h{z) 



0. 



(4.11) 



and the normalizability condition is 

dze- 3A ^h 2 (z) 



dz h (z) < oo. 



(4.12) 
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1/k 



(a) < v < 1 




1/k 



Vn 



1/k 



z 

(b) V = 1 



(c) 1< i/ < v/5/2 



(d) i/ > ^5/2 



Figure 2: Behaviour of Vh{z) for different values of v. Here, Vq = V^(l/A;). For 



£/ie radion, Vf(z) has the same behaviour with the exception that Fig. 2(c) applies 
for all v > 1 . 



In the case of study, it is only possible to obtain an analytic expression for the 
potential in the y-frame, where it reads 

2 

3e -2ky (l-fY [5v 4 k 2 (y - y s f - lQv 2 k{y - y s ) - 2z/ 2 + 5] 

VMv]) = ! ^y^r • 

(4.13) 

It is however possible to invert numerically the coordinate change (13.91) . and so 
to plot fl4.13p . Its behaviour for different values of v is shown in Fig. [2j One can 
distinguish three possible situations Q 

• v < 1 [Fig. 2(a)| In this case z extends to infinity where Vh — > 0. The mass 
spectrum is continuous from m = 0, leading to unparticles without a mass 



3 Similar potentials were considered in Ref. [24]. 
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gap. However, conformal symmetry is broken due to the occurrence of the 
scale y s . 

• v — 1 [Fig. 2(b)| z also extends to infinity but V/, — > (9/4)p 2 . This leads to 
unparticles with a mass gap m g = (3/2)p. 

• v > 1 [Figs. 2(c) and 2(d)| z s is finite and thus the mass spectrum is discrete. 
The potential diverges at z s changing sign at v 2 = 5/2, but this does not 
have observable consequences in the mass spectrum as we will see. 

Equations (14. 6 p and (14.91) do not have analytic solutions. However, for v > 1 
one can find approximations for the wavefunction in the regions near the brane 
and near the singularity. Let us first consider the region near the brane (ky ~ 0). 
Assuming ky s ^> 1 the potential (I4.13P is approximated as 

V4U - * (4.14) 

where the coordinate change is given by (I3.10p . which is approximated for v > 1 



as 



kz ~ e ky . (4.15) 

One can see that (I4.14p corresponds to an AdS metric. With this approximated 
potential, Eq. (14.91) is solved by 



h(z)\ z ~ Zo = c 1 VkzJ 2 (mz) + c 2 v kzY 2 (mz). (4-16) 

The two coefficients can be determined by the normalization and the boundary 
condition (14.111) at Zq, which yields 

c 2 _ Mm/k) fm^^ (417) 



Ci Yi(m/k) V k 

since we expect the first mass modes to be of order m ~ (z s — zq)^ 1 , and in our 
approximation k(z s — z ) ^> 1. 

Let us now move on to consider the region next to the singularity (y ~ y s ). In 
this case the potential is approximated by 

Vh \ ~ ^- 2 ^ £ 3 ( 5 - 2 ^) 1 (4 18) 

ys ~ 4z/ 4 [k{y 8 - y)} 2 " 2 /" 2 ~ 4(1 - u 2 ) 2 (z s - z) 2 ' K ' 

where we have used that, for v > 1, the coordinate change (I3.9P is approximated 
by 

P{z s - z) = -^-r[k(y s - y)] 1 ^ 2 . (4.19) 
u z — 1 
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With this approximation, Eq. (14.91) yields the solution 



h(z) = cjVmAzJ a (mAz) + c Y V mAzY a (mAz) , (4.20) 

where 

a = W^T) (421) 

and Az = z s — z. The two integration constants can be obtained by imposing the 
boundary condition at the singularity and normalizability and by matching this 
solution to the solution for the intermediate region between the brane and the 
singularity. Near the singularity (l4.2Up behaves like 

h(z) ~ c { l\Az) 3 ^ 2 -V + c f {Az) (^-m2^-2) + c W(A^)( 2 - 2 - 5 )/( 2 - 2 - 2 ), (4.22) 

where numerical factors are being absorbed in the constants Cj. We have included 
the next to leading order in the expansion of J a as we need it for computing the 
boundary condition, which reads 

e" 3A / 2 (it(z) + ^Ahj ~ c'PiAz^+W-V + 4 1} (Az)°. (4.23) 

Again numerical factors have been absorbed in d { . Note that the boundary con- 
dition is only satisfied when cy = 0, and that this condition also ensures that the 
solution (I4.22p is normalizable when v 2 < 2. 

The boundary conditions provide the quantization of the mass eigenstates for 
v > 1. In order to compute the mass spectrum for the graviton one should match 
the solutions at the ends of the space with a solution for the intermediate region. 
Unfortunately, for the parameter range we are interested in we do not have good 
analytic control for this region. However we can extract a generic property of 
the spectrum by looking at the potential Eq. (14.131) and using the form of the 
coordinate transformation Eq (I3.10p to deduce that, assuming e kys ^> 1, the 
potential has the formal 

V h {z)=f?v h {pz), (4.24) 

where Vh is some dimensionless function of the dimensionless variable pz. In other 
words we have eliminated the two scales k,y s in favour of the single scale p given 
in Eq. (13. 6p . The spectrum is therefore of the form 

m n (u, k, y s ) = p n (v) p(u, k, y s ) , (4.25) 

where the pure numbers p n only depend on the parameter v but not on the 
parameters k or y s . 



3 This behaviour can actually be seen in the limiting cases Eq. (|4.14[) and Eq. (I4.18|) . 
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1 1.2 1.4 1.6 1.8 2 

v 



Figure 3: Mass modes for the graviton, computed for ky s > 4 ^. The massless 
{n = 0) and the first 5 massive modes {n = 1, . . . , 5) are shown. 

Moreover one can find an expression for the spacing of the mass eigenstates 
by approximating the potential as an infinite well, which is valid for m 2 ^> Vh- 
The result of this approximation is 

P IX 7T 

A '"" r(i-iM — ■ < 4 - 26 > 

Note that the mass spectrum is linear (m n ~ n), and that as one approaches 
v = 1 

lim Am = 0, (4.27) 

recovering the expected continuous spectrum at this value (for v < 1 the spectrum 
is continuous too, since (14.261) is only valid for v > 1). The numerical result for 
the mass eigenvalues is shown in Fig. [3] where these behaviours can be observed. 
Some profiles for the graviton computed numerically using the equation of motion 
(14. 6 p and the boundary conditions (14.71) are shown in Fig. HI 



9 Numcrically one finds that the scaling property in Eq. (|4.25| ceases to be valid for ky s < 3, 
as discrepancies from this behavior become greater than 1%. 
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0.2 0.4 0.6 0.8 1 

z/z s 

Figure 4: KK graviton profiles in the z frame for ky s = 30 and v = 3/2, using the 
normalization f dzh 2 = 1. The massless mode (n = 0) is peaked near the brane. 
The two first massive modes (n = 1,2) are also shown. The zero mode becomes 
more peaked near the brane in comparison to the massive modes as ky s increases. 

4.2 The Radion-Scalar system 

Now we consider the spin-zero fluctuations of the system. This is 

(f>(x,y) = (f>(y) + ip(x,y), (4.28) 
ds 2 = e- 2A{v) - 2F ^ v) r ]llv dx tl dx v + (1 + G(x, y)fdy 2 . (4.29) 

With an appropiate gauge choice, the equations of motion for the y-dependent 
part of the KK modes form a coupled system with only one degree of freedom. 
The derivation of the equations is given with detail in [23], and the result is 

6" 6" 

F" - 2A'F' - AA"F - 2--F' + AA'--F = -m 2 e 2A F, (4.30) 

(/,'tp — F' — 2A'F, (4.31) 
G = 2F. (4.32) 

The boundary equations on the brane depend on the brane tension A(0). The 
precise form of the dependence can be found in [23]. At the singularity, similarly 
to the graviton case, one gets the boundary equation 

e- 4A{v) ^'(y)\ ys = 0, (4.33) 
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and the normalizability condition 

[ VS dye~ 2A M<p 2 (y) = H dzf(z) < oo, (4.34) 

J0 Jzq 

where the field has been rescaled by tp(z) = e~ 3A / 2 ip(z). 

It is convenient, as for the graviton, to use conformally flat coordinates. Rescal- 
ing the field by F(z) = e~ 3A( - z ^ 2 F(z)/4>(z), ( 14.32}) can be written as the Schroedin- 
ger equation 

- F(z) + V F (z)F(z) = m 2 F(z), (4.35) 

where 

V F (z) = lA 2 + 5 -A-At-i + 2\^) . (4.36) 



(4.37) 



The relation between the rescaled field F and the scalar field if is 

cp(z) = e 3A ' 2 

In the 2/-frame, Eq. (14.361) is given by 

e -2ky (i-A^ [3v*k 2 (y - y s ) 2 + (-6u 2 + 8v*)k(y - y s ) + 6u 2 + 3] 

V F {z[y\) = * ^ — 17 - 2 . 

4z/ 4 (y - y s y 

(4.38) 

This potential has similar form to the graviton potential, and the three situations 
presented above also apply for the radion (with the same mass gap for v = 1). A 
difference is that this potential does not change the sign of divergence for v > 1 
but, as said before, this does not have any observable consequences. 

Let us now proceed to find the approximation for the wavefunction near the 
UV brane. Taking ky ~ and ky s ^> 1 and using (14. 15j) . (14.381) is given by 

=§i, (4.39) 

and hence the solution to (14.351) is 

F{z)\ z ^ zo = dVteMmz) + CzVteY^mz). (4.40) 

The coefficients q are to be determined using the boundary conditions at the 
brane [23]. As an example, using the condition ip(y = 0) = (which we will 
use for the numerical computation) yields 

c 2 J 2 (m/k) 



ci Y 2 {m/k) 



0. (4.41) 



10 This condition holds for brane potentials satisfying d 2 X/d(f) 2 ^S> 1 [23]. 
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Figure 5: Mass modes for the radion, computed for values of ky s > 4 ™. The first 
6 massive modes (n = 0, . . . , 5) are shown. 



Next to the singularity, using (I4.19P the potential is approximated by 



Qv 2 



F\y—Vs 



fr 



Qv 1 



4z/ 4 [k{y s - y)} 2 ~ 2 /» 2 4(1 - u 2 )'< 



that gives the solution 



with 



F(z) = cjV mAzJ a (mAz) + cy V mAzY a (mAz) , 

2 + u 2 



a 



2u 2 -2 



(4.42) 

(4.43) 
(4.44) 



The behaviour of this solution near the singularity is 

F(Z) ~ c f^ z )^ 2 +m^-t) + c ^) (A ^(6^-3)/(2^-2) + c a)( A ^)-3/(2, 2 -2) > (4 45) 



Using (14.37p we can compute the behaviour of the field and apply the normaliz- 
ability condition (14.341) . 

<p{z) ~ cf\Azf/^ 2 ~V + c^Az)-^ 1 )/^ 2 " 2 ), (4.46) 

and the boundary condition f)4.33p . 

e- 3A ^^(z) ~ c f (Az)^ 2 ^/^ 2 - 1 ) + ^(Az)^^/^" 1 ). (4.47) 
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-2 L . . . , . . " . . , . . . , . . 

0.2 0.4 0.6 0.8 1 

z/z s 

Figure 6: KK profiles of the rescaled scalar fluctuations (p(z) for ky s = 30 and 
v = 3/2, using the normalization J dzip 2 = 1. The three first massive modes 
(n = 0,1,2) are shown. 

Again, the condition cy = is sufficient to ensure both the fulfillment of the 
boundary conditions and the normalizability. The scaling of the mass eigenvalues 
Eq (14.251) and the approximation (14.261) for the spacing of the mass modes also 
holds for the radion. The numerically obtained values for the masses are shown 
in Fig. [51 In comparison to the graviton mass modes of Fig. [3], note that the first 
mode for the radion is lighter than the first massive mode of the graviton. This 
can be understood recalling that the radion does not have a zero mode 0. Some 
profiles of the scalar fluctuations of the field (p are shown in Fig. EJ 

5 Other soft walls with a hierarchy 

The particular form of W, Eq. (13. ip . guarantees full analytic control over our 
solution but may seem a little ad hoc. It is natural to ask what are the essential 
ingredients of our stabilization mechanism and whether it is possible to generalize 
it to other potentials or superpotentials. The location of the singularity, and hence 

11 One can in fact show that for v — > oo, which we can only take if we resolve the singularity, 
the lightest mode tends to be massless, corresponding to the radion profile in [25]. 
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the size of the extra dimension, is given by 



Here and in the following we will assume that W is a monotonically increasing 
function of 0, i.e. W'(<p) > 0. The integral is finite whenever W diverges faster 
than W ~ 2 . 

However, the inverse volume y' 1 is, in general, not the 4D KK scale nor the 
mass gap as there might be a strong AdS warping near the UV brane. The KK 
scale is given by the inverse conformal volume z~ l (when it is finite), calculated 

as 

1^ e A ^dy. (5.2) 



It is easy to warp the geometry near the brane without affecting y s by adding a 
positive constant of 0{k) to the superpotential, leading to 

A(y) -f A(y) + ky , for W -> W + k . (5.3) 

Notice that A(y) is a monotonically increasing function of y, such that 

kz s > e kys . (5.4) 

One sees that the KK scale is warped down with respect to the compactification 
scale, a phenomenon well known in RS models with two branes [1]. In order to 
obtain, e.g., the TeV from the Planck scale we need 

*- = L iff* 2 "' (5 ' 5) 

This is not hard to achieve in a natural manner. In our model, Eq. (13. II) . it works 
so well because the exponential behavior that was introduced for large values of 
is also valid at 0(1) negative values and dominates the integral, leading to 
Eq. Q£5). 

Moreover there are many cases where z s is infinite, even though y s is finite. 
There can still be mass gaps or even a discrete spectrum, but z s is clearly inad- 
equate to characterize the energy levels. One such example is the case W = ke^ 
that leads to a mass gap. Let us be slightly more general and consider the class 
of superpotentials given by 

W{(j>) = ke^ty-faf, (5.6) 

with 0! < 0o- This superpotential is monotonically increasing for (3 > and has 
infinite z s for (3 < | , so we will assume < j3 < | . The volume y s is approximately 

ky s ~ e-*° , (5.7) 
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so, again, ky s is (mildly) exponentially enhanced when |0 O | = 0{1 — 10), 4> < 0. 
In order to estimate the spectrum, we need the asymptotic behavior of the warp 
factor in conformally flat coordinates. For large z, it is given by 

A(z)~(pz)^P, (5.8) 
where p = 0{y~ l ). The coordinate change is given by 

z ( y ) = [ e Ai - y) dy. (5.9) 



Using our trick of adding warping while keeping y s unchanged, Eq. (15.31) . we see 
that near the singularity 

z{y) - z w {y) ~ z(y)e ku ' . (5.10) 

On the other hand, adding the warping leaves A(y) nearly unchanged near the 
singularity (adding a constant ky s to infinity makes no difference). Demanding 
thus A w (y) ~ A(y) near y = y s leads to 

[My)] 1/(1 - 2/3) = [ Pw z w {y)] l/{1 - W) = M)***} 1 " 1 -™, (s.n) 

and hence 

e -ky a 

p w = pe~ ky °~ . (5.12) 

Vs 

Combining this with Eq. (15.71) we find a strong suppression of p w /k resulting just 
from 0(1) numbers. The quantity p w sets the scale for the KK spectrum in this 
case. In fact metrics of the form Eq. (15.81) have been studied in Ref. [13]. A WKB 
approximation shows that the spectrum can be approximated by 

m n ~p w n w . (5.13) 

We see that p w indeed sets the scale of the 4D masses, which are hence parametri- 
cally suppressed with respect to k. The complete superpotential that accomplishes 
a hierarchy and leads to the spectrum Eq (15 . 1 3[) is 

W{(j>) = fc(l + e^[0 - . (5.14) 

In particular the case (3 = 1/4 generates the linear Regge trajectory spectrum 
m n — Pw n appropriate for AdS/QCD models as in Ref. [9]. In this case one 
would obtain the linear confinement behaviour of e.g. p-mesons by considering an 
additional piece in our action J d 5 x^/—g e~^C nicsons . The fact that asymptotically 
A(z) ~ 4>(z) ~ z 2 guarantees that the resonances of the vector mesons follow the 
same linear law as the ones for the scalars and tensors. 

Let us conclude this section by noting that there are certainly other ways 
to obtain the mild hierarchy ky s , including moderate fine-tunings of parameters. 
What is completely generic, though, is the fact that adding warping as in Eq. (15. 3p 
leaves ky s manifestly unchanged but suppresses the masses by an additional warp 
factor e kys . 
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6 Conclusions and Outlook 



In this paper we have studied the stabilization of soft walls, i.e. 4+1 dimen- 
sional geometries with 4D Poincare invariance, that are only bounded by a single 
three-brane but that nevertheless exhibit a finite volume for the extra dimensional 
coordinate. The second brane is typically replaced by a naked curvature singular- 
ity at a finite proper distance. In particular we have studied how these soft walls 
arise in models with a single scalar field, and classified the type of models that 
can be realized as full solutions to the Einstein equations without destabilizing 
contributions at the singularities. We have proven that all admissible solutions 
result in a fine-tuning of the cosmo logical constant. 

Our main objective has been to show how to stabilize the position of the 
singularity at parametrically large values compared to the 5D Planck length. We 
have employed the superpotential method of [19, 20] and proposed a family of 
models that accomplishes this goal. Our stabilizing superpotential allows three 
types of 4D spectra: continuous, continuous with a mass gap, and discrete. The 
4D mass scale p, controlling the mass gap in the continuous case and the spacing 
in the discrete one, depends in a double exponential manner on the value of the 
scalar field at the brane, Eq (13.71) . and can thus be naturally suppressed with 
respect to the 5D mass scale k without fine-tuning. 

Next we have studied in detail the spectra resulting from fluctuations around 
our family of solutions. As with the case of the background, we have paid close 
attention to the boundary conditions at the singularity, projecting out solutions 
that would give contributions at the dynamically generated boundary. We have 
given analytical forms of the wavefunctions near the brane and the singularity, as 
well as numerical values for the lowest lying excitations and their profiles. 

We have also given a constructive recipe of how to obtain superpotentials that 
accomplish the stabilization of the hierarchy, a desired spectrum, and an "end-of- 
the-world singularity" that is consistent with the equations of motion. In a first 
step, one chooses the asymptotic (i.e. large 0) behavior of W . This will determine 
the asymptotic form of the spectrum. To ensure consistency of the equations of 
motion, the divergence should be milder than W ~ e 2 *. The different asymptotic 
forms of W and corresponding spectra are summarized in Tab. d In a second 
step, one completes W for smaller values of in such a way as to accomplish 
a mild hierarchy of the proper distance y s with respect to the fundamental 5D 
scale k, given by the simple relation Eq. (15. 5p . Notice that many of the interesting 
spectra require some kind of exponential behavior at large 0, such that this region 
does not contribute at all to ky s . At the same time one minimizes the effective 4D 
potential, Eq. (12.231) . at the brane at y = to find the vacuum value of 0o = 0(0). 

12 A classification of superpotentials giving rise to confining backgrounds was performed in 
Ref. [10]. 
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W{<j>) 


<0 2 


>0 2 
< 


e 4> 


e ^/3 
< /? < \ 


> e*05 
< e 2 * 


> ^ 


Vs 


OO 


finite 


Z s 


oo 


finite 


mass 
spectrum 


continuous 


continuous 
w/ mass gap 


discrete 






consistent 
solution 


yes 


no 



Table 1: Spectra resulting from different asymptotic forms of the superpotential. 
In the first row we give the asymptotic behavior of W{(j)), with the strength of 
the divergence increasing from left to right (> means "diverges faster than", etc). 
Second and third row show the finiteness of y s and z s , with the behavior changing 
at W ~ 2 and W ~ e^02 respectively. The third row shows the spectrum, while 
in the last one we indicate the consistency of the solution. 

In a third and final step, one adds a constant of O(k) to the superpotential. This 
adds strong warping near the UV brane, but has no effect whatsoever on the 
determination of ky s and 0o- We have shown that in this way one can warp down 
the parameter setting the overall scale for the spectrum by a factor e kys , leading 
to the desired hierarchy. 

There are a number of phenomenological applications which are outside the 
scope of the present paper but which are worth of future investigations. For the 
range of the parameter 1 < v < 2 these applications are common with two brane 
models, as RSI, but with some peculiarities. In particular graviton (and radion) 
KK modes are at the TeV scale and they can be produced and decay at LHC 
by their interaction with matter ~ h^T^", so they are expected to be produced 
through gluon annihilation [26]. Since there is no IR brane, for soft-wall models to 
solve the gauge hierarchy problem the Higgs boson (either a scalar doublet or the 
fifth component of a gauge field in a gauge-Higgs unified model) has to propagate 
in the bulk and it has to be localized near the singularity for its mass to feel 
the warping. On the other hand fermions with sizable Yukawa couplings (third 
generation fermions) should be localized near the singularity as well while first 
and second generation fermions (or near the) UV brane. As we 

have seen that the first graviton KK mode is localized near the singularity, once 
produced it is expected to decay into either Higgs or ti pairs. For v = 1 the mass 
spectrum of fields propagating in the bulk is a continuum above an 0(TeV) mass 
gap. This continuum (endowed with a given conformal dimension) can interact 
with SM fields propagating in the UV brane as operators of a CFT, where the 
conformal invariance is explicitly broken at a scale given by the mass gap, and 
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can model and describe the unparticle phenomenology. In particular the Higgs 
embedded into such 5D background can describe the unHiggs theory of Ref. [27] 
in the presence of a mass gap. In all those cases the strength of electroweak 
constraints should be an issue. Finally, the case where the spectrum is described 
by linear Regge trajectories {y = 1, (3 = 1/4) can give rise to a phenomenological 
description of AdS/QCD, similar to that of Ref. [9], where the QCD scale can be 
naturally stabilized by the scalar field. 
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